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Summary 

The highlights of the study performed under Grant NAG5-338 have been the 
conceptual design of a free-fall ground-based testing approach for gravity 
gradiometers , and the joint study and development with IFSI/CNR of the single- 
axis sensor that will bo tested at the NASA-MSFC's 300' drop-tower facility. 
Through the sponsorship of the NASA Grant, the project has reached now a point 
that we can undertake hardware development for these tests. SAO will submit to 
USAF-AEGL, on or before July 29, 1985, a proposal to continue under DoD support 
the R&D activity that started with a two-year support from NASA-GSFC. The plan 
envisions a joint testing activity NASA/USAF on the gradiometer , with AFGL 
supporting SAO for the required hardware development, and MSEC making available 
the drop-tower facility and providing assistance in personnel and equipment for 
the tests, at no cost to SAO or to AFGL. 

SAO has published three Semiannual Reports {as well as this Final Report) , 
to illustrate the accomplishments of the project; the 1st is dated December 
1983 (Fulignl and Gross! , 1983) ; the 2nd, July 1984 (Fulignl and Grossi, 

1984a) ; and the 3rd, December 1984 (Lorenzini and Grossi, 1984b). These three 
reports , together with the present one {dated July 1985) , give a fairly complete 
account of the results obtained in the course of the study. The present report 
focuses on the activity from January 1985 to June 1985. This last semester of 
the 2-year project has aimed at identifying technical approaches to inc reas e the 
sensitivity of a non-cryogenic gradiometer toward the goal of 10" 3 EU/\/Hz , 
with solutions that have the pot ential of achieving an even more ambitious 
threshold, such as 10‘ 4 EU/ \/Hz™. This goal can be achieved with a gradiometer 
design in which the proof masses are each suspended from two small arms, the 
torsion of which is directly related to the displacement of the sensing element. 
A negative-spring action, aimed at reducing the resonance frequency, is provided 
in this design by means of an external electrostatic field. This configuration 
of the instrument is also suitable for use in a tensorial arrangement. 
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1,0 INTRODUCTION 

The investigation carried out by SAO in the two-year span, 1983-to-19Q5 , 
had the primary objective of ascertaining feasibility and practicality of gradi- 
ometer schemes characterized by a sensitivity of the order of 10" J -to-10‘ 3 
EU/\/hz , that could be designed, constructed and tested in a relatively short 
time- frame, such as before the decade of the '80s is over. First, a key to the 
success of this undertaking was judged by SAO to be the direct transfer to the 
gradiometer design of high-sensitivity displacement sensing techniques, already 
developed, or under development, by the community of gravity radiation antenna 
designers (Gross! , 1981). Devices of particular interest were found to be the 
condenser probe and the microwave cavity, of the cryogenic and of the non- 
cryogenic variety. Second, it was considered crucial to the success, the abil- 
ity to test the sensitive gradiometer on the ground, in an environment free from 
man-made vibrations and seismic oscillations. The best way (perhaps the only 
possible way) to achieve this goal is, in SAO opinion, the adoption of a free- 
fall approach, using possibly (in order to save money) an existing facility, 
such as the drop-tower available at NASA-MSFC, Huntsville, Alabama. This ap- 
proach was introduced in SAO Semiannual Report #2, July - 1984 (Section 2.9, pages 
43-51), An air-tight, sealed, cylindrical container, with inside pressure 
< 10‘ 3 torr, housing in its interior the gravity gradiometer under test, would 
fall Inside an elevator cabin in rapid descent (free-fall drop is 293.8', and 
fall time is 4.3 seconds, during which data can be collected), thus simulating 
for the instrument a vibration- free , ln-vacuo, orbit-like, condition. It would 
be impractical to use bulky cryostats in free-fall tower tests, such as the ones 
that can be conducted at NASA-MSFC. This was one of the reasons why SAO, with 
NASA concurrence, concentrated its efforts on the design and construction of a 
gradiometer that uses a non-cryogenic , condenser-probe displacement-sensor, in- 
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strument. Tho plan for the construction of the overall test package is illus- 
trated in Section 2.3 (page 12-10) of Semiannual Report j|3 (December, 1984), 
The plan calls for a joint effort of SAO and IFSI-CNR, Frascati, Italy, The 
latter would construct the sensing element, using their condenser probe consist- 
ing of a rectangular plate attached to a rigid frame by crunk-shaped suspen- 
sions, all machined from a single block. The former would assemble in a suita- 
ble container the overall test package, by adding to the sensing element the 
required ancillary instrumentation such as a power amplifier to feed an existing 
modulator/transmitter for wireless transfer of data, and a gyro package to meas- 
ure tho package’s angles and angular rates during the free-fall. Data from the 
gyros would be time-correlated to the data from the gradiometer. Correction of 
errors due to attitude variations would be performed in the post-test analysis 
of the observations. 

In the last six months of project performance, SAO has done the following: 
a) further analysis of the non-cryogonic gradiometer configuration, to ascertain 
the possibility of improving the instrument sensitivity of about an order of 
magnitude; b) Investigation of the sensitivity improvement that could be ex- 
pected from an operation of the sensor at cryogenic temperatures; c) study of a 
tensorial configuration; d) evaluation of the response of a single-axis, off- 
line, gradiometer to the field of a spherical mass, in a situation such as the 
one to be encountered in the free-fall, drop-tower tests planned at the facility 
available in Huntsville, Alabama. 
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2.0 TECHNICAL DISCUSSION 

2.1 General 

The design philosophy that evolved in the course of our study favored a 
simple, non-cryogenic instrument, that could be designed and constructed in a 
two-year time frame and could bo tested on the ground, using an existing free- 
fall, drop-tower facility. Under the circumstances, it was of concern to verify 
chat, notwithstanding the instrument simplicity, the design is advanced enough 
to provide attractive sensitivity values, and is characterized by a good poten- 
tial of further improving on the threshold, without major rework of the configu- 
ration. Our most recent effort has addressed this specific issue. We found 
confirmation of the advisability to proceed in the development of a simple, non- 
cryogenic instrument that is sensitive enough to satisfy NASA-OSSA spaceborne 
gravity gradient measurement requirements, in terms of minimum detectable grav- 
ity anomaly, and at the same time is able of making use of orbital flight 
opportunities that are far more numerous than the ones suitable for bulkier and 
heavier cryogenic apparatus. for instance, the instrument of our conception 
could be flown in the second mission of the Tethered Satellite System (TSS) , 
scheduled to be launched in late 1989, while the tethered subsatellite could not 
easily accommodate the large cryostat needed with a cryogenic gradiometer. 

In the following sections, we discuss the rationale for our design choices 
and we further define our proposed free-fall testing approach. If for the 
latter all goes as expected, this approach may represent a real breakthrough in 
the development of a gradiometer for orbital flight, possibly saving years of 
efforts in the elimination of vibrations and seismic oscillations, that would 
otherwise totally mask the signals in ground-based testing. 
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2.2 Non Cryogenic Gravity Gradiometer 


Preliminary tests conducted at IFSI (Bordoni ot al . , 1985) on a modal sen- 
sor allow us to draw soma conclusions about tho performances and potentialities 
of the non cryogenic approach to gravity gradiometry. This sonnor, combined 
properly with the other identical sensors, will eventually make it possible to 
build a full tensorial instrument, The general concept of this force transducer 


has been already described in Semiannual 

Report |f3. 

Tho 

main 

featuro 

is 

that 

the proof mass is suspended through two 

small arms, 

tho 

torsion 

of 

which 

is 

directly related to the displacement of 

the sensing 

element 

and 

then 

to 

tho 


modulation of the pick off capacitors. 

Length and cross section of the suspension arms, together with tho moment 
of inertia of the proof mass with rospect to the rotation axis, determines the 
frequency of the fundamental mode of the sensor. Other modes of vibration, 
which are present, have much higher proper frequencies, and this is one of tho 
main merits of the solution that we have adopted, One further point is that 
torsion allows, fairly easily, to reach low resonance frequencies. 

This is a very important matter as far as the mechanical behavior of the 
sensor is concerned. In fact this frequency has direct influence on sensitiv- 
ity, high sensitivity requiring low resonance frequency. On the other hand, 
possible adverse dynamic conditions of the environment where measurements are to 
be performed call for a high dynamic range and strength of the instrument. 

To achieve this it would be desirable to adopt the conflicting solution of 
a high resonance frequency, or, better, to change the strength of the equivalent 
harmonic oscillator, depending on the dynamic and environmental conditions. 
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Within reasonable limits this is possible by using on artificial negative 
spring, as it is obtained by on electrostatic field. In our cose this is 
accomplished by adding two additional capacitor plntos on opposite sides of the 
proof mass, electrically insulated from those of the pick off capacitor. In 
this way half of the surface of tho sensor is used to detect the signal, and 
half to decrease tho mechanical frequency through a dc voltage regulated by a 
proper feedback loop. 

It is easy to estimate the effectiveness of the negative spring. 

As we use two plates on the opposite sides of the moving one, tho electri- 
cal forces, if of proper value, will not change the equilibrium position of the 
proof mass, and the effect of the applied electrostatic fields will be that of 
decreasing the stiffness of the mechanical spring. Such negative elastic con- 
stant is obtained by taking the raecond derivative of the electrostatic energy in 
the capacitor 

- £;(§ ov«) = oE„> (1) 

where c is the total capacitance, v the applied DC voltage, and Eo the electric 
field between plates. So the resulting spring constant, k 7 , is 

k 7 — kn — k 0 (2) 

where k a is the mechanical elastic constant. 

Equation (2) can be written as 

= «. a ( i -Po) (3) 

\ mw a ‘ J 


where m is the mass of the equivalent harmonic oscillator. 
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In thia way, In principle, the resonance frequency of tho proof macro con bo 
mada as low as one wonts, by propor choico of fl Q , 

After this brief premise, we are now in the position to ovoluato tho prac- 
tical limits of tho sensitivity of the non cryogenic instrument, in a realistic, 
even if a bit idealized, situation. 


For a gradiometer having a resulting resonance frequency of tho sensor w 7 , 
as given by equation (3), much greater than the signal frequency w„, we have for 
the sensitivity of the generic gravity tensor component P 


r J B [?~ + 2T n 

mL 2 w, [2 ^ 


( 4 ) 


providing proper mechanical and electrical matchings. This appears to be possi- 
ble in our case. This expression for sensitivity assumes that the only noise 
sources are the Internal ones, that is tho thermal source (first term in paren- 
thesis) and amplifier source (second term In parenthesis) , 


As high Q materials, like Aluminum 6061, are available, the brownian motion 
term in equation (4) can be kept sufficiently low to disregard it with respect 
to T n , the noise temperature of the electronic amplifier. So, given the mechan- 
ical configuration and parameters, tho true limiting factor is the noise of the 
amplifier . 


Wo now ask the question: given the best (from tho point of view of noise) 
available room temperature amplifier, let's find a minimum limit for sensitivity 
assuming a reasonable and practical mechanical configuration of the sensor. 

By equation (3) we can in principle decrease the resonance frequency till 
zero. This however set an upper limit to the mass of the sensor 
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m < 



(S) 


Eo must not exceed 10 7 volto/m, in order to prevent discharge, and a reasonable 
figure for c is around 1000 pE . As it was dlscuaaed previously, it is wise not 
to have a too low mechanical rosonanco frequency. 30 Hz appears a reasonable 
compromise that meets the requirements of mechanical strength, dynamic range, 
and reliable construction. With those parameters, the mass cannot be greater 
than 2.8 kg (by equation 5). 


We may now assume to be ablo to decrease at will the frequency of the 
oscillator and so to be in a position opposite to that considered in equation 
(4) , i . e . , ui| . 


In this assumption of quasi- free motion of the proof mass, the equation (4) 
for the sensitivity is modified as follows. 


~ Af ~ + 2KT n 
mL* Q 


( 6 ) 


By assuming the signal frequency at u>„ oi 0.1 rad see’ 1 , and by adopting a limit 
mass, as given by equation (5) , of 2.8 kg and a baseline of 1 m, we obtain 


r = 3 x 10- 4 EU/V^z 


(7) 


provided Q > 10 4 , which can be obtained even working at room temperature. T n 
was taken equal to 0.1 K, as it is allowed by a EET amplifier. 


The sensitivity (equation 7) corresponds to a limiting situation which is 
rather difficult to reach in practice. For example, a field of 10’ volts/m, as 
used to estimate the mass limit of the sensor, is rather unsafe, as discharges 
can be started. We are however confident that sensitivities an order of magni- 
tude worse than stated in equation 7 could in principle be obtained with more 
realistic room temperature instruments. 
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To bo specific, lot us discuss what can be obtained by pushing aa Car aa 
possiblo the porformanco of the instrument now under development at IFSI/CNR. 
The first prototype of the sensor has a mass of 0,3 kg and a resonator* frequency 
of 59 Hz, Tests on this sonsor to assoss the performance of the oloctrostatic 
field as a negative spring have shown that with an applied voltage of 50 volts, 
the frequency was decreased to 46 Hz. This result makes us confident that with 
a socond prototype, now under construction, having a mechanical resonance at 30 
Hz, an approximately 0 frequency will bo reached, 

The first gradiometor prototype will have a baseline of 0,5 m. From these 
parameters ono could oxpect a theoretical sensitivity of 

r «= 2 x 10* J EU/\/Hz (8) 

To roach effectively this sensitivity one must however solve a number of exporl- 
v,n!i ul problems, whose solutions are not at all obvious. One problem, to be 
carefully investigated, is the evaluation of noise which will bo eventually 
produced by the feedback systam. This is needed not only to adjust the scale 
factors of the sensors and regulate their dynamic range, but also to restore tho 
equilibrium position of these quasi free proof masses, We feel that it will be 
strongly reduced by the configuration of the system, where forces are exerted on 
opposite direction of tho sensor plate, providing a zero net effect, 

Another source of noise, which was not taken into account in the previous 
discussion, is the voltage and phase noise of the a.c. supply to tho bridge of 
the pick up circuit. 

Previous experience with gravity wave-antenna transducers suggests to us 
that this noise too should be strongly reduced by the geometrical configuration 
we are using. 
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2.3 Cryogenic Instruments 

In the previous section we have discussed the performance to bo expected by 
a non cryogenic instrument and we have found that the most important limitation 
to sensitivity comes out from the noise temperature of the preamplifier. 
Present day technology does not allow, at the frequencies of some kHz of inter- 
est to us, noise temperatures much lowor than 0.1 k. 

This is not true for cryogenic temperatures whore SQUIDS, or cooled GAs 
very high frequency amplifiers, in connection with microwave cavity resonators, 
con bo used. 

One further point i?t that when the amplifier noise temperature is suffi- 
ciently low, the thermal noise contribution may become important, thus requiring 
increase of 2 or decrease of temperature of the bath or both. This in turn can 
be accomplished by going to a cryogenic environment with the instrument. 

One question which deserves some consideration is that of exploring the 
possibility of using the proposed non cryogenic instrument at liquid He tempera- 
tures , with some minor modifications. 

The mechanical assembly does not pose particular problems as geometries 
similar te that under discussion have been already successfully used in connec- 
tion with a cryogenic gravity wave antenna. 

The only substantial changes required concern the pick up electric circuit 
where the FET amplifier has been replaced by a SQUID. This in principle implies 
only to change the coupling of transducer output to amplifier input. This can 
be done with a superconducting transformer. The construction of this trans- 
former, however, may be difficult in practice because a very stringent electri- 
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cal matching condition are to bo mot. In addition, in presence of high mochani- 
cal Q’s, instabilities may develop in the pick up circuit. 

In seems to us, in any case, that this line of future development is worth 
of deeper Investigations. 


2.4 The Tensorial Configuration 

The sensor we have discussod is the building block of the gravity gradiome- 
ter . A or.e axis configuration is presently under construction at IFSI/CNR. It 
is expected to have a resonance frequency of approximately 30 Hz and a mass of 
about 0.4 kg. The total length is 0,5 m. The geometry chosen is such that the 
sensor measures the off line component Tij of the gravity gradient tensor. 

The reason for this choice for the first prototype is twofold: 

(1) It has been proposed (Gross! , 1984) to make such a test in free-fall, 
in order to overcome the great difficulties of insulating the gradiometer from 
acoustic and seismic noise when testing its performance. The drop tower at 
NASA/MSFC will be used. Both at the start and at the end of the fall the 
gradiometer will experience large accelerations. Their effect on the sensors 
will be greatly reduced, and perhaps, will become unimportant, if the resulting 
inertial forces are constrained in a direction normal to that of the sensitive 
axis of the gradiometer. This, in connection with possible rotational motions 
of the instrument which can be started at the beginning of the free fall, have 
suggested to us, at least in the very first test, the above geometry. In 
addition, if the gradiometer axis remains vertical during the fall, the instru- 
ment will be insensitive to the earth field gradient. The expected response of 
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this geometry to a spherical mass simulating a gravity anomaly is discussed in 
Section 2.5 here below, 

(2) Precise geometrical relationships (both alignments and orthogonalities) 
between sensor sensitive axes and baseline directions are of primary importance. 
The geometry we have decided to use as a first step allows us to roach great 
precision in a relatively simple way and at low cost. The full one-axis instru- 
ment is in fact obtained by machining a single plato of Aluminum. It is possi- 
ble in this way to keep the centers of the sensing masses within a few microns 
with respect to the geometrical axis of the gradiometor. 

From this one-axis instrument it is easy to build a partial tensorial in- 
strument by mounting three of them along three mutually perpendicular direc- 
tions, With proper orientation of the sensor axes all the three off-line compo- 
nents can be measurod. 

To get an in-line one-axis instrument does not appear, in principle, to be 
particularly difficult. It can still be obtained from a single block of mate- 
rial by machining, but it will be more expensive. We expect also that precision 
in alignments will be a bit worse than the above one for the off-line device. 

We do not expect temperature changes to affect seriously the instrument 
performance. This requires essentially some care in attaching the gradiometor 
to its frame in order to avoid differential expansion or contraction which could 
result in deformations. Temperature gradients may be more important. Due to 
the geometry of the gradiometer, we expect that this temperature gradient will 
be possibly larger along the axis of the baseline. As the instrument is com- 
pletely made of the same material except for the insulator of the capacitors, 
the only effect caused by a strong difference in temperature between the two 
sensors 1« a variation in the value of the capacitances of the same fractional 
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order as the fractional variation of thicknoss of the insulators. This is 
estimated to be 1 part in 10* for a difference of 1°C. This may affect the 
scale factors and may require correction. 

We think however that this particular point will bo bettor evaluated, by 
performing measurements on the instrument itself and by testing its behavior as 
a function of temperature and of temperature gradients. 


2.5 Response Of One-Axis, Off-Line, Gradiometer, To The Field Of A Spherical 
Mass 


We set the mass generating the field to be tested at a certain height of 
the drop tower, and we take it as origin of the coordinates of the reference 
frame, We assume that the y axis, downward oriented, is along the vertical. 
The xz plane is the horizontal one, containing the center of mass. The gradiom- 
eter axiB is oriented parallel to the vertical and lies on the yz plane, which 
contains also the sensitive axes of the accelerometers (see Figure la) . 

In this way the yz component with respect to our system of coordinates is 
measured. It is related at every point during the fall to the components T 
along the principal axis by a rotation around the x axis. This gives for the 
output signal, S, of the gradiometer 

S sa ^Tyy “ * sin 20 Pyx cos 20 


and , being r y * = 0 , 


S 



20 


or, by eliminating 0 , 
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S „ 2CM ^ 1 ,1 , r 

Ro 3 R„ Ll+ (h/R 0 ) 2 J 

This is zero when the gradiomoter crosses the xz plane and has, in absolute 
value, two maxima at h - ±R 0 /2 , corresponding to 

S m ± 0,57 GM/R 0 3 . 

If R 0 = 2m, M - 10 kg, S corresponds to 4 x 10' 3 EU , 

The form of the signal, us a function of time, is shown in Figure 2. 

In the test it is desirable to contain the signal as seen from the gradiom- 
eter in a narrow range of frequencies starting from zero. We can obtain an 
estimate of this range by considering that the duration At of the signal is of 
the order of the time the gradiometer takes to span a distance R„, as it is 
clear from Figure lb, 

At =? R 0 / 2gho 

where h 0 is the height above the xz plane at which the free fall is started. We 
must then keep R 0 as large as possible and h e small. 

With a mass of 10 3 kg we get a signal of 10" 1 EU at 7 m. By assuming 
h„ = 10 m, we get approximately for the frequency within which the bulk of the 
signal Is contained 

Af = 2 H* . 

An enhancement, due to the form of the signal, is to be expected at a 
frequency f 0 = l/(2t D ) , where t„ is the time needed to travel the distance 
between the two peaks, which gives f, a 1 Hz. 
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Figures 3, 4 and 5 give the spectrum of the signal at the output of the 
gradiomoter. The spectrum was obtained by FFT'ing the simulated instruments ' s 
output in the geometric situation of Figure 1. Figure 3 is the most informative 
of the three Figures and provides the Power Spectral Density (PSD) of the sig- 
nal. It can be seen that the peak of PSD is at about 0.5 Hz. Figure 4 and 
Figure 5 give the real and the imaginary part of the signal spectrum. 

Appendix I contains the software code that was used in computing the graphs 
of Figures 3, 4 and 5, Appendix II contains the print-out of the calculations. 
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Figure 1 a 

Geometry of free-fall measurements. 
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Figure 1 b 

Graph of Signal S as a function of distance h 
(See Figure 1 a^for meaning of symbols R Q and h) 


The analytical expression for S is: 


S = 


2 G M 

* 


h 

% 


m 


1 

1 + (h/Ro) 1 
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Figure 2. Signal at the output of the gradiometer as a function of time. 










3,0 CONCLUSIONS AND RECOMMENDATIONS 


We conclude on tho advisability of a follow-on program based on the con- 
struction by SAO of a non-cryogonic gradiomatbr package that usas the displace- 
ment sensor developed at IESI/CNR, and on the froo-fall testing of such a pack- 
age at the drop-tower facility available at NASA/MSEC, Huntsville, Alabama. It 
is our understanding that NASA-OSSA and USAE-AEGL have already formulated a 
joint activity plan, with AEGL providing financial support to SAO for the in- 
strument hardware effort, and NASA-OSSA providing support for the in-house ac- 
tivity at MSEC. We recommend that such a plan bo implemented without delay; we 
actually expect that this follow-on phase will start in October 1905. If this 
is the case, there will be time enough to prepare for a 1909 flight, or ovon for 
an earlier flight. 

It appears that NASA Grant NAG5-338 has made it possible for SAO to achieve 
enough progress in gradiometer design, that hardware development is now on the 
verge to be initiated. If the events in Eall 1905 will confirm this expecta- 
tion, tho Grant would have fully accomplished its scope. 
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Appendix I 

Software Code for the Computation of the Spectrum of the 
Expected Signals in Eree-tfall Experiments 
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GRADSPEC.FORJ2 17*34 PAGE 1 

PROGRAM GRAOSPEC 
IMPLICIT REAb*B(A-H f n-Z) 

C GRADIOMETER RESPONSE FUNCTION , 

C COMPUTE THE FFT AMPLITUDE. 

C PRINT THE RESULT, WITH THE FRF.OURNCY, 

C NOTE THIS T5 THE SQUARE ROOT OF THE PSD, 

C WE SUBTRACT OFF THE 11, C, COMPONENT BEFORE FFT* TNG, 

C IN > ANAL t ' THERE HAY HE PROBLEMS DUE TO THE HANNING WHEN WE 

C IGNORE THE D.C, POST-FFT AS being COMPONENT FFT(t). 

PARAMETER UMAX « 2048 , IHAX2 * 1 025) t IMAX2«IMAX/2+l 

DIMENSION ST(IHAX),5T0C(IMAX),P0W{IHAX2) , 

1 FREOCIHAX2) 

PARAMETER (ISIZE»1024) IHALF THE MAXIMUM N 
COMPLEX* I 6 Z ( ISIZE+ l ) 

COMMON /FFTCOM/Z I HAKE FULL FFT ACCESSIBLE OUTSIDE 

LOGICAti HANNING (FLAG IN COMMON - SET Til DO HANNING, 

COMMON /FFTCOM2/HANNING 

PRINT *, ' ( PI* T**2 - P2)/P,Y * 

' PRINT *, * ENTER Pt,P2,PJl' 

READ *, PI , P2, P3 

PRINT *, 1 TIME INTERVAL! THIN , TMAX 1 

READ *, THIN,TMAX 

PRINT », 1 NUMBER OF INTERVALSl * 

PRINT *, • (MUST BE EVEN * i**K IS MOST EFFICIENT)* 

READ *, NPTS 

H ■ (THAX-TMIN)/(NPTS-t) 

DO I a l, NPTS 

T * THIN 4- ( I* I ) *fl 
X » (Pi*T*T-P2)/P3 
ST ( I ) a X / C (t,0+X**2)**2,5 ) 

END DO 

WRITE( 12,101) Pl,P2,PJ,TMIN f THAX,NPTS 
101 FORMATC PARAMETERS I*, IF10,3,» INTERVAL! 1 , 2F10.1,* I POINTS*, 

l 15) 

C PRINT *, 'ENTER NUMBER OF POINTS IN THTER VAN ' 

C READ *, NINT 

NI, NT a NPTS 

PRINT *, ' ENTER 1 IF HANNING IS TO BE DONE' 

READ *, IMAN 
IF (IHAN.EO.l) THEN 

HANNING a .TRUE. 

WHTTE( 12 » * ) > HANNING DONE* 

else 

HANNING = .FALSE. 

WRITEC 12, * ) ' HANNING *NOT* DONE* 

ENDIF 

PRINT ENTER 1 IF D.C, LEVEL SHOULD BE SUBTRACTED! ' 
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GRAUSPEC,F0RJ2 29-MAR-P95 17(31 

READ *, I DC 
IF (inr.EO.l) THEN 

CAM, 5URnc(ST,NINT,STflC) 

CAM, FFTAHP(STDC,NtNT,POW) 
WRTTF.M2,*) ' D » C , LEVEL SMITH ACTED 1 

ELSE 

CAM, FFTAMP(ST,NINT,POW) 

WHITEC12,*) * D.C. LEVEL SIIDTRACTED' 

EUDIF 

NP a NINT/2*! 


EQUAL INTERVALS 


gf^ a l,/( NI»T«HJ^tjAS5l)HE EC 

~FRFQUENCI INTERVAL!'# DF# 


WRITEC12,*) 

DO I b 1 i HP 

freq ( i ) ■ < i- 1 j *nr 

F.ND no 


(HERTZ) ' 


215 


HRITEC12,*) ' F,P, REAL, COMP!,' 

WR HE( 12 # 2 1 5) 

1 (F«EO(I),POWm#Zm, 1 = 1, NP) 

FORHATCO', // 

1 (F15,7,5X,G16,8,5X,2G16,0) ) 


END 


SUBROUTINE SURDCC A#N#8) 
IMPLICIT REAL*fi(A-H,0-Z) 
dimension A(i),nn) 

AHEAN a 0. 
on 1 a 1 , N 

AHEAN a AHEAN + A(T) 
END DO 

AHEAN a AMEAN / N 

DO I o 1,N 

BCD = A(I)-AHEAN 
END DO 


RETURN 

END 
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FFTAMP,F0R|J2 2<1-MAR-IRH , > 1 7 t 3S PACK t 

SUBROUTINE FFTAMP(A,N,XAMp) 

IMPLICIT HF,A!.*A(A-H,n-7,) 

OILERS ION A(I),XARP(1) 

PARAMETER (ISIZKb|024) JHALF THF MAXIMUM N 

complex*^ zcisizkm) 

COMMON /rr TCOH/Z JHAKC FULL FFT ACCESSIBLE DUTS10F 

LOGICAL HAN NT MO I FI.AG IN COHHON - SET TQ DO HAN NT NO, 

RATA HANNING/, TRUE,/ t DEFAULT NET ON 

COMMON /FFTCOH? /HAN KING 
INTEGER |WK(6*ir>tZE+150) 

REAU*H WK(6*I5TZF.+ t5R) ,AA(2*ISIZE) 

PARAMETER (PT33,l<U5R2f.S!5BR1 t n,PI2 = 2*Pl) 

INPUTl 

A( 1 , ,N) SEQUENCE TO BE TRANSFORMED 

N SIZE OF SEQUENCE — 

MUST BE EVEN f. POSITIVE, tip TO 1024, 

OUTPUT! 

XAHp(l,, N/2+1) AMPLITUDE OF RESULTING DFT 
XPUASKt " ) PHASE " * " 

COMPUTES A PROPERLY NORMALIZED D,F,T, DF THE REAL INPUT A(). 

OUTPUTS, RATHER THAN REAL AND COMPLEX COMPONENTS , THE 
AMPLITUDE! NOTE THAT PHASE IS LOST. 

BEFORE COMPUTING DFT, APPLY A HANNING 

FUNCTION MULTIPLIER TO REDUCE "LEAKAGE" ( SEE BRIGHAM, "THE 
EAST FOURIER TRANSFORM", P.140)! THIS REDUCES FAKE BASELINE AND 
SHARPENS UP LINES AT FREQUENCIES "BETWEEN" THE DISCRETE 
FREQUENCIES OUTPUT, BUT SLIGHTLY BROADENS LINES PRECISELY AT AN 
OUTPUT FREQUENCY. TO GET POWER IN A LINE, SUM THE AMPLITUDES 
ACROSS THE LINE BEFORE SQUARING, 

USE IMSL ROUTINE FFTRCt ) ! LINK WITH SYS6 IMSL J IMSL9P/LIB . 

N IS NOT REQUIRED TO BE A POWER OF 2, THOUGH MORE EFFICIENT IF SO, 

C APPLY HANNING! 

SCALE ■ PI2/N 
DO I b 1 , N 

IF (HANNING) THEN 

AACI) = A(I)*0,*5*( 1-COSC (I-1)*SCALE I ) 

ELSE 

AAU) » A ( I ) 

ENOIF 
END DO 

C TRANSFORM! 

CALL FFTRC(AA,N,Z,1WK,WK) 

C GET AMPLITUDE! 

INP ■ N/2 ♦ l 
SCALE * 2*n/H 
DO I o 1 , INP 

XAMP(i) b ARS(ZCt))*SCALF. 

F.ND DO 


RETURN 

END 
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Numerical Results Obtained with the Software Code 
Illustrated in Appendix I 



Pago 31 


OKiG&iiL FAC*.-" 

OR FOOii CU^Vh'V 


PARAMETERS J 5.000 10,000 7, 

HANNING On HE 
n,C, tiEVF.E SIJOTOACTEn 
FRFOURNCY TNTEOVAltl 0 . 2495 1 t 7 J H7500000 
F,P, PEAL, C^RPli 

000 IMTFRVAhl 
{HERTZ) 

0,000 4,000 

utinormalized 

f 

m 

Real [F] 

2-m 

0,0000000 

0.40614849E-0I 

10.402571 

n.oootiofiooF+no 

0,2495117 

0.6640953 IE-01 

-16,070007 

-5.50/7591 

0,4990214 

0 ■ 7 19557 R2E-01 

17,607370 

-6,971561 9 

0 , 7 <05352 

0.647O6461E-01 

-5.678R5R6 

15,5/9164 

0*9900469 

0 , 4609909 1 E-0 1 

-5,5049707 

-10.669736 

1,2475596 

0, 31343031E-01 

7,7074557 

2.2 10R651 

1 , 497070 J 

0,1066761 0 E-0 1 

-4,1245741 

2,5704453 

1 ,7466020 

0.1 17H0176F.-01 

0,67176413 

-2,9400066 

1,9960918 

0.692O1520E-02 

0,97730122 

1.4776005 

2.2456065 

0 , 393200 I 2E-02 

-0,9920016b 

-0.17076590 

2.4951172 

0,21 9698 92E-02 

0,45471407 

-0.31056406 

2,7446299 

0,122620 4 3E-0? 

-0.40106965E-01 

0. 31 1 35623 

2,9941406 

0 • 05676602E-0 1 

-0,10650490 

-0,1 .1009621 

3.24J6523 

0 , 1582 1648K-03 

0. 9 122 3(14 1 E-0 1 

0,93741 M'JE-02 

3,4931641 

0, 19460239E-03 

-0,36741 34HE-01 

0, 316441 J.lK-01 

3,7426750 

0 , 97 40 64 92E-Q4 

0 , 1 (1000 1 60F-02 

-0,2491 2279 ( '-0 1 

3.902JR75 

0 , 5600 40 79E-04 

0.94057570E-02 

0, 1 075OR 4 RE-0 1 

4.2416092 

0,2R192607E-l)4 

-0,7lH0S6B6E-02 

0,7272960 1 F-03 

4.4912J09 

0 , 1 10 1 n 1 97H-04 

0 ,25904 l 96F.-0? 

-0 , 209 34 06 1 E- 02 

4,7407227 

0.96227R72E-05 

0.0265R010E-04 

0 , 24620497E-02 

4.9902344 

0.30653551E-05 

-0,73 190459E-03 

-0,21770b 15E-0 3 

5.2397461 

O.241H01 I2E-05 

0,515770 15E-03 

0,1426 1797.K-01 

5,4992570 

0.27I742R4E-05 

-0,16711 270E-03 

0,642506 JlE-,03 

5.7307695 

0,62653 10BE-06 

-0. 1 6925U54E-04 

0, 1594969HE-03 

5. 90 B 2.0 1 3 

0, 1 1944467, K-05 

0.557027 12E-04 

0, 100661 R t E-03 

6.2177910 

0.92940080E-06 

-0.14771 651 E-04 

0,235 37 205E-03 

6.4973047 

0.74072486E-06 

0,1 01 J557IF.-04 

0,| R9J5449E-01 

6,7369164 

0,74724 156E-06 

0 , 1 900246 1 E-05 

0,1 9120491 E-03 

6.9063291 

0,6790240 IE-06 

-0.3060341 4E-05 

0,1611 RR74E-03 

7.235R308 

0 ,5733466 IE-06 

0 . 2, 2 6 9 3 R 4 JE« 0 5 

0,1467591 9 E-0 3 

7,4051516 

0, 524O0O97E-06 

-0,53610215 E-06 

0,1341431 0E-03 

7.734R613 

0 , 469768 1 BK-06 

-0,1 490707 6E-06 

0,1201 3 266E-0 3 

7,9941750 

0 , 4792060 1 E-06 

0,27780671 E-06 

0, 109H°6«7E-03 

8.2310067 

0,391 237 72E-06 

-0.12577144E-06 

0. tnO15670E-OJ 

0.4831904 

0.35737390 E-06 

0 ,437257 3HE-07 

0,914077 2 9 E-04 

0,7329102 

0,32005401 E-06 

0,2461 190RE-07 

0,01901 8 22E-04 

9 , 992421 9 

0.301 3754TE-06 

-0,7096nOH1E-00 

0,771521 20E-04 

9.2119316 

0.27765416E-06 

0,17 176 190F.-07 

0 ,7 1 07°464Fi-04 

9.4914453 

0.25617675E-06 

0 ,64292 1 70E-00 

0, 65632449 FJ - 04 

9.7309570 

0.23716731E-06 

0.62 185576E-0H 

0.607148 16E-04 

9.99O46R0 

0,21 90501 7.E-04 

0,76211 6O6E-O0 

0,66201677 F',-04 

10.2299005 

0.20418377E-06 

0,54201 440E-08 

0.5227103JF-04 

10.4794922 

0, 1R99fii54F,-06 

0 , 646000 9 3F, -00 

0,406 30 1 54F’.-04 

10,7290019 

0, 1 7 70 3224E-06 

0,4971001 IE-00 

0 , 45 320352E-O4 

10.9706J 46 

0, 165247 32E-06 

0 , 440 3909 1 E-OfJ 

0,4230331 IE-04 

11.2200273 

0,154486 ), 0E-0 6 

0,41 161269E-O0 

0,195481 H5E-0 4 

11,4775391 

0.I4463677E-O6 

0, 37 4 0075 9E-0H 

0, 3702701 2F -04 

U. 7270500 

0,1 156070IE-06 

0.34U 1126 E- OH 

0, 14715 395E-04 

11.9766625 

0,177 1I7.B4K-06 

0,31 63209 JE-OB 

0, 17,59200 7 F >04 

12,2260742 

0,1196007 4E-06 

0, 29O7R7OOE-O0 

0, 3O63026HE-O4 

12,4755859 

0,11264595 E-06 

0,26041 539E-0H 

O.70H37364E-O4 

12,72.50977 

0.1 Oft 151 71 E-06 

0,240011 qgE-OH 

0.771 747 0®F— 04 

1 2,9746094 

0. 1O014564E-06 

0 , 2294 7 463 E-OB 

0. 256.17204E-04 


points 512 
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I 1,2241211 
I3.473632H 
1 2 . 77 J 1 44*« 
13,9726563 
14,2221680 
14,4716797 
14,7211914 
14,9707011 
IS, 2202148 
15,4697266 
15,7192381 
15,9687500 
16,2182617 
16,4677714 
16,7172052 
16,9667969 
17,2163066 
17,4650203 
17,7151320 
17,9648418 
IB, 2141555 
18,4638672 
10.7133789 
18,9628906 
t 9, 2 ( 24023 
19,4619141 
19,7114258 
19,9609375 
20,2104492 
20,4599609 
20,7094727 
20,9584844 
21,2084961 
21 , 458O07B 
21,7075195 
21,9570313 
22,2066410 
22,4560547 
22,7055664 
22,9550781 
23,2045898 
23,4541016 
23,7036(31 
23,9531250 
24,2026367 
74,4521484 
24,7016602 
24,9511719 
25.2006836 
25,4501953 
25,6997070 
25,9492188 
26,1987305 
26,4482422 
26,6977539 
26,9472656 
27,1967773 
27,4462891 
27,6958008 
27,9451125 
28,1948242 


0,940H4064E-I)7 
0,8 9 4261 8 91XJ7 
0,04636039007 
0,00181417007 
0,76033449007 
0,7216619 3007 
0,68556293007 
0,651 H26BHO07 
0,62026351007 
0 ,59070064007 
0,56290220007 
0 , 5 3696652007 
0,51 752475007 
0,4895396011-07 
0,4679(1409007 
0.44752051007 
0, 4282996761-07 
0,410159550(17 
0 , 39302535007 
0,37682820007 
0,361504 B 400 7 
0,34699708007 
0,311251 30007 
0,32021 00800 7 
0.307H5176E-07 
0,2961 1018007 
0,28495432007 
0.27434H0SE-07 
0,26425707007 
0,25465271007 
0.24550376E-07 
0,23 67 8415 007 
0, 228*6895(^07 
0,220534921X17 
0,21 296042O07 
0,20572526007 
0 , 1 98 B1 059007 
0,1 921 908OE-O7 
0,18587345007 
0,1798191 4007 
0,17402144007 
0,1684660511-07 
0,16314272007 
0,1580 37171X17 
0.15313903E-07 
0,14843704007 
0,14 192.174007 
0, 13958748007 
0, 13542032007 
0,1314 1 40611-07 
0.12756O96E-07 
0.123H5372O07 
0. 12028546007 
0.11684967007 
0,11364022007 
0,11035131 1X57 
0.10727745E-07 
0,10411 345E-Q7 
0,101 45441 E-07 
0,98695675000 
0,960328 35E-08 


0.21270008O0H 
0,19739401 E-OH 
0,1834161 4008 
0,17071 3H4E-O0 
0,1 5906254E-08 
0. 14H39074E-0H 
0,1 3859902OOH 
0 , 12060 | 64000 
0,12132315000 
0, 1 1369527K-0H 
0, H)6657B8i;-0B 
0,1001560 JE-00 
0.0414 3841E-09 
O,fl0h75fe70E-O9 
n,BJ41 IJ77O0V 
0.78b2|951E-09 
0, 7416994?, E-09 
0 , 7OD20B65E-O9 
0,6617 I046E-09 
0,625792 34E-09 
0,5922b J49E-09 
0.56O95280E-09 
0,531 68 637E-09 
0 , 504 J06461.-09 
0,4786688 2 E-09 
0,451641 84E-09 
0,432 1O656E-09 
0.41095224E-09 
0, j9l 07920E-09 
0,37 2 39523 E-09 
0,154816231X19 
O.3J026434K-O9 
0, 12266H89E-09 
0,30796 J8 5E-09 
0 , 29408944E-09 
O f 28P9B967E-09 
0,26861 329E-09 
0 ,2569 1 306E-09 
0.245H4507E-09 
0,2 15 3h87 2E-09 
0 , 225446R9E-09 
0,2 16043791', -09 
0,2 0712041 E-09 
0,1 9867006E-09 
0,1 906 4 1 1 6E-09 
0,18301 548K-09 
0, 1 767696 2K-0 9 
0.16B8H050E-09 
0, 16232761F--09 
0. 1660914fll>09 
0,1601 5382E-09 
0 , l 4449770E-09 
0,139 1 07 37E-09 
0,1 3396782E-09 
0, 12906542E-09 
0. 124JH705E-09 
0,1 199207 1 E-09 
0,11 56545JK-09 
0,1 1157B60E-09 
0 . 107682 39E-09 
0 , 1 0 J95624E-09 


0,2421 3520E-04 
0.22891104E-04 
0.21666H76F.-04 
0,205264431X14 
0, 194645631X14 
0,1 84745451X14 
0, 1755041 11,-04 
0 , 1 ()6 0676H1'-O4 
0,158 1874bl>01 
0,151219 16E-04 
(l , 1 44 1 2 J441X14 
0,1 J746J431X14 
0.1 31206 341.-04 
0, J25J2214R-04 
0, 1 1978345K-04 
0, 1 145652bE-04 
0,109644701.-04 
0, 105UOU84K-O4 
0,1 00b 1 449E-04 
0,964680191,-05 
0 , 92545239E-05 
0,88831 252E-05 
0,Bb3123331>05 
0,81 975H27E-05 
0,7081 0051 !>05 
0,75B0420bt>US 
0.729483051X16 
0,702311 OOE-05 
0,676500151X15 
0.65191099E-05 
0, 6284H96 3E-05 
0,606167 42E-05 
0, 584H0O5UE-O5 
0 , 564569 39E-05 
0,5451 7B67E-05 
0,52665bhbE-05 
0 , 50 8965 1 01X15 
0.49202B94E-05 
0, 47583b04i:-05 
0, 460J3690E-0b 
0 , 4 45404H BE-0 5 
0,431275151X15 
0, 4 I 7 64537K-05 
0 ,4045751 5E-05 
0,192035921X15 
O,3BOOOO07E-O5 
0.36844479E-05 
0 , 157 J4394E-05 
0 , 34667602E-05 
0, 33642OO0E-05 
0, 32665606E-05 
0,3 1 7(J6bb2K-05 
0, 3079 3077E-05 
0,2991351 5K-05 
0,290b62961X)5 
0 ,282499 J5E-05 
0 , 27 4b 10271X15 
0 , 26704 24 4E-0.6 
0, 259723 101.-05 
0 , 25266 09 3K-06 
0, 245844061X15 


2d, 444)359 

0,93461 7 1 41.-00 

0,10il)V2051. 

09 

1 ), 23926199) 

0 5 

2d , o *» J R 4 7 7 

0, 909 M) 3561.-08 

0,969H,)97911 

10 

0,73290459). 

05 

2B.Q43 1504 

0,HR57K907I>0« 

0,9J 1)5)641' 

t 0 

0,22676223)1 

05 

29,1 9287 1 1 

0 , 862600761-00 

0,91)5078031: 

10 

0.720H7579). 

05 

29,4423820 

0,040187071,-00 

0,8759147711 

10 

0,21 50H66DE 

05 

29,691X945 

0,016501 5971-00 

0,847195245 

to 

0,7095 164 1). 

06 

29,94140ft) 

0 . 7975208 7)1-01! 

0,81 9bb498K 

ID 

0,204167 JOE 

05 

30, 1 9091 Bl) 

0,77723478)1-00 

0.79)258065 

10 

0, 19B9721UE 

05 

30,4404297 

(1,757691 66)1-00 

0,7679257011 

10 

0,19)94346); 

05 

30,60994)4 

0, 7 3067 3 l 51,-00 

0.74)61 JObll 

10 

0,1890747 3)1 

05 

30,93945)1 

0,7201 542OH-O0 

0 ,77.027074El 

10 

0,18435947)1 

0 5 

31,101)9640 

0,707)109 1)1-00 

<1,697(1583611 

10 

0.17979160)1 

05 

31,4)84756 

0, 68502050 K- OH 

0,6763236*11 

to 

0.17536527E 

05 

31 .6H79B03 

0.66826146K-0R 

0,655635)06 

10 

0, 17107 493)1 

05 

31.9375000 

0 , 6520 | 2901,-08 

0,6)574701B 

10 

0,166915)0)1 

05 

32,1070117 

0,6 362551 6K-0B 

0,6166242211 

10 

0, 162881)2), 

05 

32,4365234 

0,b2096942K-(tfl 

0,5982)073)1 

10 

0,15896 II I7E 

05 

32,6860352 

0,606137701, -08 

0,5805395811 

10 

0, 15517125)1 

05 

32.9355469 

0.591742B4K-0H 

0,56351 06 HE 

10 

0, 15148617 E 

0 5 

33,1050506 

0.57776B44E-0B 

0.547121401; 

10 

0, 14790872): 

05 

33,4345703 

0,5641 98B 3E-0R 

0,531 3)91 7E 

11) 

0,I444)490K 

05 

33,6040020 

0,55101902)1-08 

0.516138H4E 

10 

0,14106087)1 

05 

33,9)35936 

0.5 J82146B)1"0H 

0,5014971b); 

10 

0,1 J778296E 

05 

34.1631055 

0,52577217)1-08 

0,4873863 OR 

10 

0, 1 J45976b)l 

05 

34,4326172 

0.51367B21K-08 

0 , 47 37 R563E 

10 

0,13150162)1 

05 

34,6021209 

0. 5019203911-08 

0, 4n(>674 35E 

10 

0,128491 62E 

05 

34,9316406 

0. 49048664 K- 00 

0 , 44802944E 

to 

0,12556458 E 

05 

35,1011523 

0.479J6544K-D0 

0 , 4 J68 )45bE 

to 

0, 1227 1755E 

05 

35.4306641 

0 , 4685 457 6E-0B 

0,^2406421 E 

to 

0,1 1994771E 

05 

35,6801750 

0,458 0 1701 E-OH 

0.41271090F, 

10 

0,117252)5)1 

Ob 

35,9296075 

O.447769O5K-0B 

0,4 01749)4 11 

10 

0,U46?8HBK 

0 5 

36,1791992 

0.437792J 4E-0B 

0 , 391 1 6740R 

to 

0,1 1207479)1 

0 5 

36,4207109 

0.42007696K-00 

0.38094692E 

in 

0.10958770E 

0 5 

36,6702227 

0,41 B6 1455)1-08 

0, 37107802E 

10 

0, 10716533E 

05 

36,9277344 

0,409)963111-00 

0,161 5423 IE 

10 

0,1048054bK 

05 

37, 1772461 

0,4004139fll>08 

0, 35232754E 

10 

0.1025059UE 

05 

37.4267570 

0.3TI65962B-08 

0, 3434240611 

10 

0, 1 (1026486)1 

05 

37,6767695 

0,30.11256711-00 

0, 3348140611 

10 

0, 9B0B0 I 59E 

06 

37,9267813 

0,374B0465K-0H 

0,3264923011 

10 

0 , 9594999 1 E 

06 

38,1752930 

0,3666896811-00 

0.31H44J67E 

to 

0,9 387255HE 

06 

30.424R047 

0,3587739311-08 

0, J10b5797E 

to 

0,91846125)1 

06 

30,674)164 

0. 35 1050H 811-88 

0, JOJ12604E 

10 

0,89869025)1 

06 

38,9230201 

0,34?i5| 42811-08 

0.295H4026E 

10 

0.H7939655E 

06 

39,1733390 

0,3361580911-08 

0.28H79002E 

to 

0 ,8606647u)l 

Ob 

39,4220616 

0,32897650(1-1)8 

Q.2H196158E 

10 

0,H42179HS)1 

(1b 

39,672)6)3 

0,3219639411-08 

0 , 27535649E 

to 

0,82422768)1 

06 

39,9210750 

0,3151 ( 50011-08 

0.26H95869E 

10 

0,80669441): 

06 

40,1713867 

0,3084245111-08 

0.26276371E 

to 

0,78956675)1 

Ob 

40,4208984 

0,30) 8874711-08 

0 , 2567627 JE 

to 

0,772HJ193)1 

Ob 

40,6704102 

0,2954990611-08 

0.25095078E 

10 

0,75647759)1 

0 6 

40,9199219 

0,2892546311-08 

0, 24537048 11 

to 

0,7 40491 BbE 

06 

41,1694336 

0.28314971E-08 

0 , 2 398630H E 

10 

0.724H6324E 

06 

41,4189453 

0,277 1 7996K-0 8 

0.23457607E 

10 

0.70958U70E 

06 

41,6684570 

0,271 )41 2)11-08 

0,229451 33E 

10 

0.6946H54H 

06 

41,9179600 

0.2656294811-08 

0.2244H266E 

10 

0,68001 147)1 

Ob 

42,1674805 

0.26004084H-08 

0,7196677111 

10 

0.66570454E 

Ob 

42,4169922 

0,254571 5411-08 

0,714990)711 

10 

0,65170315)1 

Ob 

42.6665039 

0,2492179711-08 

0, 2104697 OE 

to 

0,6379980 IE 

06 

42.9160156 

0,2439766)5-08 

0 , 20607 04 HE 

10 

0,62458017)1 

06 

43,1656273 

0,2)88411 311-08 

0.2018214HE 

to 

0, 61H409b)l 

06 

43,4150391 

0,2338172011-08 

0 , 1 976B65bll- 1 <1 

0.59H57202E 

06 
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or^.iv . . i ’* 
OF irOi55i “ 


43, 6645508 

0.278H9768K-08 

0,1*3 36797011 

Ml 

0,58596525k 

do 

43,91 40675 

0,22406751 k>0H 

0,189792261. 

10 

0,57 Jnl2H Ik, 

0 6 

44,11.35742 

0,219 1 J075OOH 

0,18b0209ji; 

10 

0, 56I50719E 

06 

44,4130B59 

0.21470151 K-00 

0,1823568811 

10 

0.54964100E 

06 

♦ 4»bli25977 

0,21015906^-1)8 

0, 17HB0598E 

10 

0. 53800737k; 

06 

44,9121094 

0.205'70267U-08 

0,1153676111 

10 

0,52659HH4E 

06 

4b.16M.2U 

0,20li3179f:-()« 

0, 17201082E 

10 

0,515409 J8E 

06 

45,411132(1 

0 , 19704 J88E-08 

0,1687627411 

10 

0 ,5044 32 34k' 

06 

45,6606445 

0,1928 J652K-08 

0,1656097811 

10 

0,49 3661 50k, 

06 

45,9101563 

0, 1H0707 35E-IJU 

0,1625493 IE 

10 

0, 48 30908 1 K 

06 

46,1596680 

0, 18465407H-08 

0,159577001, 

to 

0,47271 44.IE 

06 

46,4091797 

0,18067 4 48E-OB 

0, 156693 3 1 E 

10 

0,462b2668E 

06 

46, 65669 1 4 

0, 1767664 1H-08 

0 , 1 5 3B91 40E 

10 

0, 45252205B 

06 

46,9067031 

0, 17292787K-0H 

0,151 1722111 

10 

0,442b9522K 

06 

47,1577148 

0,16yi5664i;-08 

0, 1485360611 

10 

0,43 304 IO0E 

Ob 

47,4072266 

0.165450931.-OH 

0, 1459676 5E 

10 

0,42355438k: 

06 

47,6567383 

0,161 80H78K-08 

0,1434842611 

to 

0,4t423047E 

Ob 

47,9067500 

0.158 2 2B 34E-IIH 

0, 1410654911 

10 

0,405064555' 

06 

48,1557617 

0.1.54107016-08 

0, 13H72125E 

10 

0, 39605201b 

06 

48,4052734 

0,151 2454611-08 

0,13644.16211 

Ml 

0,387188 3 8E 

Ob 

48.6547852 

0. 147839591,-08 

0, 1342305511 

10 

0, 37H46934E 

06 

48.9042969 

0.1444885411-08 

0,1320865511 

10 

0,3698V0bhE 

Ob 

49,1536066 

0,141 1 907 2fcl~0H 

0,130001 37E 

10 

0, 3bt 44824K 

Ob 

49,4033203 

0,1779445711-08 

O.12797019E 

10 

0,3531 380911 

06 

49.6528320 

0,1347405711-08 

0, 1260159111 

10 

0,3449563 Jb 

06 

49,9023438 

0,1 316012511-08 

0, 1241 1054E 

to 

0.33689920H 

06 

50.1518555 

0 , 1 28 50 1 1 BE-OH 

0,1 22260 1 HE 

30 

0, 32896301 b 

06 

50.4013672 

0.1254469511-08 

0, 1204669411 

10 

0,321 1441 mC 

06 

50,6508780 

0. 12241721E-08 

0,118725236 

10 

0,3134392611 

06 

50,9003906 

0,1194706711-08 

0, U703534E 

10 

0, 30584482b’ 

06 

51,1499021 

0, 1165459 IE-08 

0,1 154000 3E 

10 

0.298J5757E 

06 

51,3991141 

0,11 36618711-08 

0,1 1 38071211 

10 

0 , 29091 429E 

Ob 

51,6489258 

0,110817 1211-06 

0,1122674911 

to 

0,28 3691 H.lk' 

llh 

51.8984375 

0,10801060^-08 

0, 1 10773 19E 

10 

0,2765071 Jt 

06 

52,1479492 

0,1052410911-08 

0, 1093279811 

10 

0 , 2694 1 7 1 9E 

06 

52,3974609 

0, 1025074711-08 

0,1079?J2)K 

10 

0,2624191 IK 

06 

52.6469777 

0 . 9980860611-00 

0, 1065625JE 

to 

0 , 2555 1003b 

06 

52,8964844 

0,971 4 3426E-Q9 

0,1052455611 

10 

0, 24868717E 

06 

51,1459961 

0,9451086211-09 

0,10 197092k' 

10 

n,?41947HlE 

Ob 

53,3955078 

0,9190987911-09 

0,1027365811 

10 

0,7 3528929K 

06 

53,6450195 

0, 0971945911-09 

0,1015406011 

to 

0,228709O2E 

06 

53,894531 3 

0 , 067986 1 2E-09 

0, 10038620E 

to 

0.22220445E 

06 

54,1440430 

0,8 4 2B6366E-09 

0,9926824711 

tt 

0,2157731 (IE 

0 6 

54.3935547 

0.01801771 67-09 

0,9818805211 

It 

0,20941 253K 

06 

54,6430664 

0 , 79 7 4 J897F1-09 

0.97144531E 

tt 

(I. 20312038b 

06 

54,8925781 

0.7691 JB31K-09 

0,961 3531511 

11 

0, l 96H9429E 

06 

55,1420898 

0, 74604681)11-09 

0,9516758511 

tl 

0,1907 3190H 

06 

55,3916016 

0,7212156911-09 

0,9422424711 

1 1 

0,10463122k: 

Ob 

55,6411173 

0,6076163911-09 

0,9 3 320822E 

It 

0,1785898 (IE 

0 6 

55.8906250 

0,6742404811-09 

0 ,9245046 1 E 

It 

0,172605566 

06 

56,1401367 

0,6510796811-09 

0,9161488511 

1 1 

0,16667640k,' 

Ob 

56,3896484 

0, 6281 2507E-O9 

0.90806B22E 

tl 

0,160800228 

06 

56.6391602 

0,6053710611-09 

0.90027677E 

tt 

0, 15497490E 

Ob 

56,0886719 

0,5828074211-09 

0, B92889B IE 

11 

0, 14919870K 

Oh 

57,1381836 

0,5604272411-09 

0,8857345511 

tl 

0, t43469J7kl 

(16 

57.3876953 

■9 , 5 3822291 E-tlO 

0,8788897211 

n 

0,137 7 8S07E 

Ob 

57,6372070 

0,5161069811-09 

0.B7234346E 

tt 

0, t 32143U7K 

«b 

57,8867188 

0,4943120811-09 

0,0660960011 

Lt 

0.12654J89E 

Oh 

58.1762305 

0.4725909611-09 

0.86012799E 

tt 

0,t20983 29E 

06 

58.3857472 

0,451016491,-09 

0,8544521711 

tl 

0, U546022E 

Ob 

58.6352539 

0,4295816111-09 

0.84902205E 

It 

0,109 ** 7289 H 

Ob 
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vi^uau l :,iS. : f v 

OF POOR QUALITY 


50 , PB47656 
59,1342773 
59,3037991 
59,6333000 
59,9020175 
60,1323242 
60.3010359 
60,6313477 
60,0000594 
61,1303711 
6t , 3790020 
61, 6 7. 91945 
61.0769063 
62.1784100 
62.3779297 
67.6274414 
62,6769531 
63,12646^0 
63.3759766 
63,6254003 
63,8750000 


0.40H77936R-09 
0. 307 IO2nBF><>9 
0, J6604539K-09 

Q,345inoi*f>09 

0.32426059G-09 
0,30 352Q09K-09 
0,26267717 (■> 09 
0,26211 OJ9K-09 
0,241«2637f>09 
0.22141976K-09 
O.2O1O7032K-O9 
0. 10O7q777K-O9 

o, 16057 

0.14039461K-09 
0.12025970G-09 
0, 10016107E-09 
0.BO092670E-10 
0.6004B1D3E-10 
0.40022234E-10 
0.2000B1 19R-10 
0.31054152E-13 


0,84 3B6679E-1 1 
0 , 9 J0 1 1 0 1 9E» 1 1 
0, 6 J44O220E"! 1 
0,nj009664E-U 
O.026O469OE-1 1 
0,B2727719E»1 1 
0.81H70565E-11 
0 , 0 1 63 1990E* 1 1 
O.P123B731E-1 1 
O.0O957O9OE-11 
0,60701 175K-11 
0.80470977E-U 
0.R0766974K-H 
0, 60092252 K«1 1 
0, 79910849 E ”11 
0.79H07827E-11 
0.796B5702E-U 
0.7959H550E-U 
0,7950 5007E” 1 1 
0.79474205E-H 
0,79498630011 


0, 10451952006 
0,99090 J 1 ?E-07 
0, 9 370761 9(-;»07 
0, 0 0 34 564 IK-67 
«,H30107U)K”07 
0,77701 141K-07 
0.72415275K-07 
0,67151 459E-07 
0,6 1 9O0O67E”O7 
0,566h3462E-07 
0.5V476O51K-O7 
n , 462042 JOE”07 
0,411 0641 2K-07 
0, 34041020E-01 
0.307H64H2E-07 
0,2564 12 3 3K"07 
0 , 2050 3 722K*07 
0, 1 5372364 E» 07 
0.10245669K-07 
0,5 122072 3E»06 

n,ooooooooE+oo 



